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Abstract 

We study a two-parameter extension of the cosmological standard model ACDM in which cold 
dark matter interacts with a new form of dark radiation. The two parameters correspond to the 
energy density in the dark radiation fluid AA^uid and the interaction strength between dark matter 
and dark radiation. The interactions give rise to a very weak “dark matter drag” which damps 
the growth of matter density perturbations throughout radiation domination, allowing to reconcile 
the tension between predictions of large scale structure from the CMB and direct measurements 
of eg. We perform a precision fit to Planck CMB data, BAO, large scale structure, and direct 
measurements of the expansion rate of the universe today. Our model lowers the ^-squared relative 
to ACDM by about 12, corresponding to a preference for non-zero dark matter drag by more than 
3(7. Particle physics models which naturally produce a dark matter drag of the required form include 
the recently proposed non-Abelian dark matter model in which the dark radiation corresponds to 
massless dark gluons. 
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I. INTRODUCTION 


Cold dark matter (DM) has long been an important ingredient of the cosmological stan¬ 
dard model ACDM. Evidence for its gravitational effects exists on a range of length scales 
from galaxy rotation curves and gravitational lensing to large scale structure (LSS) and 
the cosmic microwave background (CMB). Particle physics models of DM generically also 
predict non-gravitational interactions of DM. Such interactions could be responsible for de¬ 
termining the DM abundance via thermal freeze-out and would alter the clumping of DM at 
small scales. Until now all attempts to observe non-gravitational interactions of DM more 
directly have only yielded upper bounds. The searches include “direct detection” of DM 
collisions with a target in the laboratory, “indirect detection” of DM annihilation in regions 
of high DM density in the spectrum of photons emanating from such regions, and “collider 
searches” which seek to produce and observe DM in particle collisions in the form of missing 
energy. 

Another possible route to learning about the nature of DM is through precision cosmolog¬ 
ical measurements. Precision fits to the CMB, large scale structure, and several probes of 
the background expansion based on baryon acoustic oscillations (BAOs), supernovae lumi¬ 
nosity, measuring the current expansion rate H 0 , etc., are sensitive to detailed properties of 
the dark matter. For example, hot dark matter in the form of active neutrinos may not con¬ 
tribute to more than 1.9% of the total dark matter density today, which turns into a bound 
on the total neutrino mass Tim v < 0.21 eV (95%CL) |T] - or Tim v < 0.12 eV when using 
also Lyman-a forest data from quasar spectra [2]. The dominant component is compatible 
with the assumption of cold dark matter (CDM), or possibly warm dark matter (WDM), 
but with strong bounds on the velocity dispersion of warm particles in order to avoid a cut¬ 
off in the matter power spectrum on scales where precise data are available [3jJ. Similarly, 
consistency of the global fit limits the density of additional particles which may contribute 
to cosmic radiation to an equivalent number of neutrino species N e g < 3.7 (95%CL), while 
their masses satisfy model-dependent bounds, of order rn x < 0.38 (T v /T x ) 3 eV (95%CL) 
for one thermal species x [T] . Recent Planck data also requires standard neutrinos to be 
free-streaming, c 2 viscosity = 0.331 ± 0.037 (68%CL) p. 

The “cosmic concordance” ACDM model is impressively successful and some of its pararn- 
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eters have now been measured to permille accuracy. However, with the increase in precision 
of the measurements, there have also been indications in the data for effects which are not 
well described within ACDM. Perhaps most significantly, there is tension between the value 
of erg predicted within ACDM with parameters fit from the CMB and BAO, and the value 
of <jg from more direct measurements of LSS with various techniques (CMB lensing, galaxy 
weak lensing, cluster mass function, etc.) This tension could be caused by systematics in 
the astrophysical data, or may point instead to new physics affecting dark matter, because 
cr 8 (roughly speaking the amplitude of matter fluctuations at scales of 8 Megaparsec) is pre¬ 
dominantly determined by the growth of fluctuations in the DM component, during matter 
and radiation domination. By taking into account neutrino masses, or by introducing extra 
massive relics, it is possible to fit CMB and BAO data with smaller erg values, but this is 
done at the expense of increasing the tension with measurements of H$ from local redshift 
data, and of degrading the agreement between the matter density fraction Vt m inferred from 
CMB data and from other cosmic probes [Ij. 

In this paper we take the “ag-problem” seriously and interpret the discrepancy as evidence 
for new physics affecting the DM. In particular, we investigate whether the tension between 
the CMB fit and the direct measurements of erg may be resolved by including a very weak 
drag force between DM and radiation, which acts during radiation domination. Such an 
interaction would dampen the growth of density perturbations in the DM fluid, and therefore 
reduce the predicted matter power spectrum. 

There have been several previous investigations of modifications to the matter power 
spectrum in the presence of new DM interactions, e.g. Refs. mna. In these setups, the 
new interactions tend to be very strong at some scales and negligible at others, leading to 
threshold features in the power spectrum - similar to WDM models. The data disfavor 
such features and put strong bounds on the interaction parameters, leaving little room for a 
reduction of the power spectrum at 8 Mpc scales. There have also been attempts to solve the 
problem with decaying dark matter (DDM) [TT1 fTTTj . In that case, the strongest constraints 
come from CMB data: in order to avoid a huge “late integrated Sachs-Wolfe” effect, the 
DM lifetime must be so large (see e.g. [13]) that a 8 can only be slightly reduced. Hence, 
current DDM studies do not find very significant evidence for a non-zero decay rate. 

One way to solve the ag-problem without introducing other tensions is to consider a drag 
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force which is both very weak and which results in a momentum transfer rate between the 
DM and radiation which scales with temperature in the same way as the Hubble expansion 
rate during radiation domination, i.e. proportional to T 2 . A model which predicts this kind 
of drag force and damping of density perturbations was recently proposed in [15| (BMS). 
With such a scaling, the drag can act equally on density perturbations which enter the 
horizon at different times during radiation domination. This leads to a smooth suppression 
of the cold dark matter power spectrum at all scales during radiation domination. Since 
the growth of matter fluctuations is primarily modified during radiation domination (rather 
than matter domination), there is no enhancement of the late integrated Sachs-Wolfe effect in 
such models, and CMB data do not provide strong bounds on the interaction parameter. In 
Section |IT[ we make this suggestion concrete and define the drag force, calculate the resulting 
momentum transfer rate, and discuss its qualitative effects on matter perturbations. 


In order to add an interaction for the DM with radiation, it would appear most economical 
to couple DM to either photons or neutrinos. However, in both cases, it is difficult to get 
the desired scaling with temperature. For example, a momentum transfer rate proportional 
to T 2 is not compatible with Compton scattering off photons which scales as T 4 , Coulomb 
scattering off electrons which scales as T 3//2 or weak interaction scattering off neutrinos which 
scales as T 6 . In addition, any new interactions of especially photons but also neutrinos are 
tightly constrained by the global fit to the CMB MUD]. Following the concrete model 
in BMS we therefore introduce a new self-interacting “dark radiation” (DR) component to 
the energy density of the universe which may be described as a perfect fluid. The drag 
force arises from scattering between DR particles and the DM. Thus our proposal is to 
generalize ACDM with two parameters. One corresponds to the energy density in the fluid 
describing the new radiation component (parameterized and normalized like an effective 
number of extra neutrinos AAfluid); and the other is the DM-DR interaction rate r 0 (this 
parameter gives the interaction rate today, but since its scaling with temperature is known, 
it characterizes the strength of the new drag force at any time). 


The precise definitions of these parameters and their effect on the CMB and matter 
power spectrum are given in Section [TTJ In Section III, we review the example proposed 
by BMS of a particle physics model with a non-Abelian dark gauge group, and summarize 
the calculation of the momentum transfer rate. We also offer an alternative model with a 
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massless dark photon coupled to the DM and massless fermions. Section IV contains the 
main results of this paper. We perform a fit to CMB, LSS and BAO data and find a strong 
preference for non-zero interactions between the DM and DR. The best-fitting models have 
To — 1.6 x 10“' Mpc -1 ~ 1.6 x 10“ 21 s -1 and AiVfl U id < 0.7, and the minimum % 2 improves 
by 11.4 over that of the minimal ACDM model with AiVfl U id = 0 and To = 0 . The fit shows 
a ~ 3.7cr preference for a non-zero value of the drag coefficient. A non-vanishing DM drag 
of the kind that we are proposing is also found to be compatible with large values of the 
Hubble rate, as measured for instance by [16] |T7] (when including such data, the minimum 
X 2 improves by 12.7). 


II. GENERALIZING ACDM WITH DARK RADIATION AND DARK MATTER 
DRAG 


A. Modified cosmological perturbation equations 


We propose to add a new component of self-interacting dark radiation to ACDM. The 
radiation is comprised of relativistic particles with a self-scattering rate which is fast com¬ 
pared to the Hubble rate during radiation domination. Such radiation can be described as 
a perfect fluid with speed of sound c 2 = 1/3 and no viscosity. 


We assume that the dark radiation was in thermal equilibrium with the Standard Model 
particles early in the evolution of the universe, and that it decoupled during freeze-out of the 
dark matter. Then the temperature of the dark radiation today will be of the same order 
as the photon temperature. The photon temperature is expected to be higher because the 
photons inherit the entropy of heavier particles in the Standard Model during their freeze- 
out. The main effect of the DR (other than providing drag to the DM) is to contribute 
to the expansion rate of the universe, because it increases its average energy density. We 
choose to parameterize the energy density in the DR in analogy to an effective equivalent 
number of neutrino species 


AiVflUid = iVdr 


/Tdr 

\Tv 


4 
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(bosonic DR), 
(fermionic DR). 
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( 1 ) 



Motivated by the specific particle physics models described in Section III we will study 


values of AiV fluid ranging from 0.07 to 1 (see section IVA). 


The dark matter in our model is comprised of particles which become non-relativistic long 
before matter-radiation equality and therefore have negligible kinetic energy density. We 
take the dark matter to have negligible self-interactions and parameterize its contribution 
to the energy budget of the universe with oj dm = fi dm h 2 as usual. 


The drag force between the DM and DR can be parametrized by the linear coefficient T 
of friction which a non-relativistic DM particle of velocity v experiences as it propagates 
through the thermal bath of radiation 

v = —aT v , (2) 


where the dot here and in the following represents a derivative with respect to conformal 
time and a is the scale factor. The coefficient T depends on the temperature of the DR. In 
any specific model it can be computed (see Section 3) from the rate of momentum transfer 
due to collisions of the DM particle as it travels though the DR. For the case of interest - 
CDM - the velocities in Eq. ([ 2 ]) are non-relativistic and satisfy c 2 v 2 T dr /M where c is 
the speed of light and M is the DM mass. 

We use the formalism of Ma and Bertschinger [18} and write the coupled evolution equa¬ 
tions for density and velocity perturbations of the DM, <5 dm and 9 dm , and the DR, 5 dl . and 
9 dr . In Conformal Newtonian Gauge, the equations for the DM and DR overdensities in 
Fourier space are 


Aim — 

— 9dm + 30 

(3) 

0dm = 

--0dm + GT(6» dr - 9 dm ) + k 2 0 
a 

(4) 

Air = 

4 

— g^dr + 40 

(5) 

9 dr = 

e s f+ei,+ 3 l Pdm a m m « dr ). 

4 4 p dr 

(6) 


Here p dm and p dr are the average energy densities of DM and DR, respectively; and 0 and 
0 are the scalar metric perturbations in Conformal Newtonian gauge. Notice the absence 
of all higher moments of the DR perturbations, they vanish for perfect fluids. 

The purpose of displaying these equations is to draw attention to the drag terms propor- 
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tional to T which couple the two velocity equations. They represent the drag forces which 
result from collisions between the particles in the two fluids. To gain a rough understanding 
of what these terms do, note that during radiation domination the coefficient in the equation 
for the dark radiation is suppressed by the small ratio pdm/pdr- Therefore the main effect 
is the drag due to the DR on the DM. Note that the clock in these equations is set by the 
Hubble rate, which scales like T 2 during radiation domination. Therefore, if we require that 
the effect of the drag term is small at any instant, but uniform over a long interval of time, 
it must also scale as T 2 . We will assume that this is the case, and show how to motivate 


this behavior with concrete particle physics models in Section III 


After matter-radiation equality, the Hubble rate decreases more slowly, proportional to 
T 3 / 2 , whereas the drag continues to be proportional to T 2 . Thus the effects of the drag 
become less important after equality. We use T 0 to denote the value of the drag coefficient 
extrapolated to today. The drag coefficient at any other temperature is then 

2 

(7) 




where T 0 = 2.7255 K is the current CMB temperature. 


B. Effects on the CMB and LSS spectrum 

We implemented the above model in the Boltzmann code CLASS 1 US ESQ] . Very few 
modifications of the public version of the code are required for this model. We implemented 
the new equations in both the Newtonian and Synchronous gauge, and checked that we get 
exactly the same results in the two gauges. The only difference is that in the Newtonian 
gauge, we can run with a density of ordinary non-interacting cold dark matter cu c dm set to 
exactly zero, while in the synchronous gauge we must set it to a negligible but non-zero 
value, e.g. cu c d m = 10~ 10 , since the latter gauge is by definition comoving with the CDM 
component. The code assumes natural units (c = 1) and expresses conformal time and 
Fourier wavenumbers in Megaparsecs (Mpc). Hence T 0 is naturally expressed in inverse 
Mpc. It can be converted to inverse second by multiplying by 0.97 x 10~ 14 Mpc/s. 

1 github.com/lesgourg/class_public or class-code.net 
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FIG. 1: Ratio of the dark matter density perturbation 5dm for an interaction rate To = 2 x 
ICR' Mpc -1 ~ 2 x ICR 21 s^ 1 over the same perturbation in the standard non-interacting limit, in 
the Newtonian gauge, as a function of conformal time, and for six representative wavenumbers. The 
interaction rate causes a suppression of 5 dm inside the Hubble radius, efficient especially during 
radiation domination, and continuing during the beginning of matter domination (the vertical 
dashed line shows the time of equality between radiation and matter). Apart from To, the two 
cosmological models share the same parameters, including AAfl u id = 0.21. 


In Section IV, we will find that models with a rate of the order of T 0 — 2 x 10“ 7 Mpc -1 ~ 
2 x ICR 21 s~ x provide the best fits to the data. Figure [l] shows the evolution of 5dm for such 
a value of To normalized to a ACDM model with To = 0. In the figure, AfVfl U id = 0.21, but 
we will later show that the effect of r 0 and MVj ui d are not strongly correlated. The figure 
shows the evolution of six different wavenumbers between k = 10~ 3 Mpc _1 to k — IMpc -1 . 
We see that the growth of DM fluctuations is suppressed roughly between r ~ 2n/k (time of 
Hubble crossing) and the beginning of matter domination (roughly until r ~ 2500 Mpc, the 
time at which the ratio of radiation to matter density is of order 0.1). Later on, i.e. deep in 
the matter dominated regime and during A domination, the growth curves are horizontal, 
showing that the growth rate is the same as in the ACDM model. 

Figure [2] shows, first, the effect on the temperature and polarisation CMB spectra of in- 
















FIG. 2: Residual of the temperature (left) and A-polarisation (right) power spectrum in several 
extended models compared to the minimal ACDM model. Two models have ordinary decou¬ 
pled cold dark matter, but either free-streaming (blue) or self-interacting (green) extra relics 
with respectively AIV e g- = 0.21 or AAg u id = 0.21. The text explains which quantities have 
been kept fixed in these comparisons. The last model (red curves) shares the same parameters 
as the latter model (green curves), excepted that the DM-DR interaction is switched on, with 
To = 2 x 10~" Mpc -1 ~ 2 x 10~ 21 s _1 . The boxes show the binned error bars of the Planck High 
Frequency Instrument 2015 data , which covers l > 30. All models are well within the error bars 
of the Low Frequency Instrument, which covers i < 30. 

creasing A N e g from zero to 0.21, in a model with extra free-streaming relics (e.g. relativistic 
sterile neutrinos), and in our model with self-interacting dark radiation AiVj^ = 0.21. For 
a useful comparison, the redshift of radiation/matter and matter/A equality are kept fixed 
by appropriately scaling (Ud m and H 0 . The baryon density Ub and reionization optical depth 
are constant. This transformation absorbs a significant part of the total effect, except for 
perturbation effects (in particular, the gravitational drag exerted on the photons by the 
DR), and an enhanced Silk damping effect. While the Silk damping effect is the same in 
the free-streaming and self-interacting models, the perturbation effects are not. In the free- 
streaming case, the extra radiation has very smooth perturbations, inducing extra damping 
of the CMB spectra. In the self-interacting case, DR features larger perturbations, which 
boost the CMB spectra through gravitational coupling. In the latter case, the Silk damping 
and DR gravitational drag effects tend to compensate each other, and the net variation of 
the CMB spectrum is smaller than in the free-streaming case. For that reason, we expect a 
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FIG. 3: Residual of the matter power spectrum P(k , z = 0) in the same extended models as in the 
previous figure, compared to the minimal ACDM model (see the caption of figure [2] for details). 

weaker bound on AiVfl UK ] in the case of self-interacting DR. In both models, the .E-mode po¬ 
larisation spectrum is also affected on large angular scales (/ < 40). Indeed, an enhancement 
of the radiation density changes by a small amount the thermal history (different freeze-out 
value of the free electron fraction at the end of recombination, and primordial Helium abun¬ 
dance inferred from BBN). In order to maintain the same reionization optical depth, the 
redshift of reionization changes slightly, and the shape of the low -l polarisation spectrum is 
affected at the level of a few percents. However, because of cosmic variance, this effect is of 
small relevance when fitting CMB data. 

In Figure [2j the last (red) curves show the additional effect of switching on the DM-DR 
interaction, with r 0 = 2 x 10“' Mpc -1 ~ 2 x 10“ 21 s“ 1 , all other parameters being fixed 
like in the model with self-coupled dark radiation and ANq^ = 0.21 (green curves). The 
impact of the DM-DR interaction is hence given by the comparison of the red and green 
curves. Overall, this impact is small, since for the same value of To, the CMB is affected at 
the level of ~ 2%, while the matter power spectrum is affected by 20 to 30% on the range 
of scales most relevant for <r 8 . This follows from the fact that in any model in which DM 
fluctuations evolve on a time scale set by the Hubble rate (rather than some shorter time 
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scale imposed by microphysics), there is an effective gravitational decoupling between DM 
and photon fluctuations (2U 22j. Models with DM-DR interactions can sometimes violate 
this condition, and generate “dark oscillations” [7] with a period T <C if -1 . However, in 
our model and for the range of parameters in which we are interested, dark oscillations 
remain negligible, as clearly shown by Figure [l] Hence the effective gravitational decoupling 
still holds in good approximation, and primary CMB fluctuations are weakly affected by 
modifications in the DM growth rate before photon decoupling. Still, they can be affected 
by modifications in the DR growth rate. When the interaction rate r 0 is switched on, DR 
clusters differently. This propagates to photons through the usual DR gravitational drag 
effect, i.e. the gravitational interaction between photons and dark radiation. In the CMB 
spectra, this type of gravitational dragging effect is effective on scales slightly smaller than 
the Hubble scale at any given time, and shifts the acoustic peaks in phase and amplitude. 
This explains most of the oscillatory features visible in the red curves of Figure |2j that start 
near l ~ 200, i.e. on scales which are slightly smaller than the Hubble rate at the time of 
photon decoupling. The red curves also encodes smaller effects depending slightly on r 0 , 
like the early ISW effect, and lensing by small scale structures. 

Figure [3] shows the effect of AA^y and r 0 on the matter power spectrum evaluated today, 
P(k,z = 0), for exactly the same models as in Figure [2} As discussed in BMS [15], scales 
crossing the Hubble radius during matter domination are unaffected by the drag effect, while 
smaller scales are reduced. The suppression depends on the time spent by each mode inside 
the Hubble radius during radiation domination. In this example, in the range of scales 
contributing to a 8 (k ~ 0.2/rMpc' 1 ), the matter spectrum is suppressed by about 20%, as 
needed for solving the cr 8 problem. 

III. EXAMPLE MODELS AND CALCULATION OF A fluid AND T 0 

In this Section, we briefly review the model of non-Abelian dark matter and dark radiation, 
more details can be found in BMS [15]. We also define an alternate model with a dark photon 
and give expressions for A#j ui( i and r 0 for both. 

In the non-Abelian DM-DR model the dark matter is a Dirac fermion which transforms 
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as the neutral component of an S77(2)-weak triplet. It’s couplings to the Standard Model 
are identical to the pure “wino” DM in some supersymmetric models, however our model 
does not have supersymmetry. The DM particle also transforms in the fundamental “iV” 
representation of an SU(N) dark gauge group. Therefore it interacts with “dark gluons” 
of the SU(N) gauge group. For phenomenological reasons we will be interested in dark 
gauge couplings of order ad ^ 1CT 8 and gauge groups of size N — 2,3,4. Since the dark 
gauge coupling is so small, the dark gauge interactions do not confine until length scales 
much larger than the size of the visible universe. Therefore the dark gluons correspond to 
weakly-interacting massless particles, they form the DR of the model. 

In [15] it was shown that the DM in this model obtains the correct abundance from 
thermal freeze-out for masses 1.2 TeV (N=2), 1.0 TeV (N=3) and 0.9 TeV (N=4). The 
model safely evades current DM direct detection bounds, it will be within reach of future 
indirect detection experiments, and its DM could be discovered at a future 100 TeV collider. 
Since the dark gluons are so weakly coupled they do not play a direct role in DM detection 
phenomenology. However, at temperatures above the DM mass, the dark gluons come to 
thermal equilibrium with the Standard Model so that they obtain the same temperature. 
After DM annihilates into Standard Model particles and freezes out the dark gluons decouple 
from the SM and evolve with their own temperature. The temperature of the DR, T^, may 
be estimated from entropy conservation as T^/T 7 = (#*/5 , * ec ) 1//3 where g * is the effective 
number of Standard Model degrees of freedom and g* ec is the number of SM degrees of 
freedom during the decoupling of the dark radiation. 

Using this relationship and assuming that the DR decoupled with temperatures in the 
10-80 GeV range, we compare the energy density in the DR fluid with the energy density in 
standard neutrinos to obtain the effective number of neutrinos [15] 

AAfl uid = 0.07(1V 2 - 1) . (8) 

An important difference between radiation in the form of neutrinos and the dark gluons is 
that the dark gluons can be described as a perfect fluid with zero viscosity. This is because 
the rate of self-interactions of gluons - even for couplings as small as ad ~ 1CD 10 - is faster 
than the Hubble rate at any time from before nucleosynthesis until today. 

The drag coefficient T is computed by considering scattering of dark gluons off the dark 
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matter. The leading logarithmically enhanced contribution comes from a t-channel Feynman 
diagram and was computed in [15] • By resumming hard thermal loops the calculation can 
be improved to also obtain finite pieces, the results are in [23.; [23] • Here, we only need the 
logarithmically enhanced contribution 


r 0 = 


today 


1 

1 

to 

1 


' —a 2 d log a d ' 


'1.2 TeV' 


3 


2.0 x 10- 16 


l 

_l 


~ 1.9 x 10 Mpc 

where in the last line we plugged in representative values for the parameters, N 
a d = 10 —8 5 , M x = 1.2 TeV, and converted to units of inverse megaparsec. 


(9) 
( 10 ) 
= 2 , 


We close this Section by giving an alternate model for DR with couplings to the DM which 
predicts values for AA'Aiid as low as 0.07. In this model the DM is a Dirac fermion wich 
couples to a massless dark photon, the gauge boson of a dark U (1) gauge group. In addition 
to the DM and the dark photon the model also includes massless fermions which are charged 
under the dark 1/(1) (in an anomaly free representation). For simplicity, we consider a single 
Dirac DM particle with charge 1 and Nj species of massless Dirac fermions with charge 
gdr- The DM is assumed to couple to the Standard Model through the weak interactions or 
through the Higgs portal. Then it can obtain its abundance from standard thermal freeze-out 
(he. it is a WIMP) with a mass of order the weak scale. The massless fermions and the dark 
photon constitute the DR of the model. The dark photon equilibrates with the Standard 
Model at temperatures above the DM mass for the range of gauge couplings of interest. 
However, whether or not the light fermions equilibrate as well depends on the precise values 
of their charges. We find that for q dv ~ 1/3 the coupling of the light fermions to the SM 
thermal bath are too weak to bring them to thermal equilibrium before the decoupling of 
the dark sector from the SM. In this case they could remain significantly colder than the 
SM until after the decoupling of the dark photons from the SM. Then they would give a 
negligible contribution to the effective number of radiation degrees of freedom. However, for 
light fermion charges greater than 1 the light fermions equilibrate before decoupling of the 
dark radiation. After DM freeze-out the interactions between the SM and the DR become 
too infrequent to maintain thermal equilibrium and the two fluids decouple. As in the non- 
Abelian model, the DR is described by a perfect fluid with zero viscosity at all temperatures 
relevant for the growth of density perturbations and the CMB even in the case of smaller 
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light fermion charges. 

We obtain 

AiV fl uid = 0.07(1 + ^) (11) 

for the case of gdr ~ 1 and 


AWiuid — 0.07 


for gdr ~ 1/3, and values of Afl ui d ranging between the two limiting cases for charges 


q<h 


< 


1/3. For the drag coefficient we obtain 


T 2 

-i 1 d 


2vr 

To = N f qi-al logaj 1 ^- 

X 


~ 1.8 x 10 7 Mpc 


today 


1 

pV/dl 


' -a/log ot d 


"1.7TeV" 


2 


2.0 x 10- 16 


l 

_1 


( 12 ) 
1 < 


(13) 

(14) 


IV. FIT TO CURRENT DATA 


A. Data and methodology 

We use the code MontePython [25j to fit the model to currently available cosmological 
data. We split the data into four categories: 

• CMB: we use the Planck 2015 TT + low-1? likelihood from Ref. [26]. 

• BAO: we use measurements of ZV/Lirag at z — 0.106 by 6dFGS [27], at z = 0.15 by 
SDSS-MGS [28] . at z = 0.32 by BOSS- LOWZ [29], and anisotropic BAO measure¬ 
ments at z = 0.57 by BOSS-CMASS-DR11 [29]. 

• LSS: we use three probes of Large Scale Structure: the Planck 2015 lensing like¬ 
lihood pH], the constraint <T 8 (fL m /0.27) a46 = 0.774 ± 0.040 (68%CL) derived from 
the weak lensing survey CFHTLenS m, and the constraint (T 8 (fl m /0.27) a3 ° = 
0.782 ± 0.010 (68%CL) from Planck SZ cluster mass function [32]. The latter con¬ 
straints should be taken with a grain of salt, because they have been inferred under 
the assumption of a ACDM model. However, our model produces a featureless matter 
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power spectrum on the scales probed by these experiments, so these constraints are 
probably valid to a good approximation. 

• H 0 : we occasionally also use the constraint Hq = 73.8 ±2.4 km/s/Mpc (68%CL) from 
Riess et al. [16J. Direct measurements of the local Hubble rate by e.g. US HZI have 
been questioned recently by the community, with the concern that systematic errors 
might have been underestimated. However we will use it only in order to show that 
our model is well compatible with such high values of the Hubble rate. 

We did several MCMC runs with various combinations of these data sets, for the 6-parameter 
ACDM model, and for our 8-parameter model (with a free effective number of dark gluons 
AiVfl U id, and a dark matter-dark gluon interaction rate r 0 , expressed in the code in inverse 
Mega-parsecs). We define the ACDM in the same way as the “base model” in Planck 2013 
and Planck 2015, including two massless and one massive neutrino species with m = 0.06 eV, 
and assuming an effective neutrino number N e g = 3.046. We keep exactly the same settings 
in the model with interacting dark matter and dark gluons; in that case, the density of 
usual CDM is set to zero, while the density of non-abclian dark matter is parametrized by 
Wdm = f^dm h . 

We impose flat priors on the parameters of our model: {(n fe , a; dm , AAfl uid , r 0 , H 0 , A s , 
n s , T reio }. Only the lower edge of the priors on AAfl uid , To an d r re i 0 are relevant. For the 
interaction rate, we just require r 0 > 0. We impose a prior r reio > 0.04 on the optical depth 
to reionisation 2 . Finally, for the dark gluon density parameter AA% lid , our non-Abelian DM- 
DR model predicts discrete values AAfl ui d = 0.21,0.56,1.05,... for N = 2,3,4,..., while the 
second model discussed in section [lH| predicts AA fluid = 0.19, 0.32, 0.44,... for Nf = 1, 2, 3,... 
in the case of large light fermion charges and AAfj uid = 0.07 for small charges. Hence we 
stick to the theoretical prior AAfl uid > 0.07 in all our MCMC runs. This means that the 
standard ACDM model (with To = 0 and AA% jid = 0) is not a special point in the parameter 
space of our “dark matter drag” models. This is not a problem from the point of view of the 
statistical analysis, as long as we provide a way to evaluate the goodness-of-ht of the “dark 
matter drag” model compared to ACDM. For that purpose, we performed some companion 

2 This is done in order to avoid the limit T re j 0 — > 0, which is unphysical given the residual ionisation 
fraction after recombination. 
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runs with ACDM and the same combinations of data; in each case, we report the difference 
between the minimum y 2 of the two models 3 . 

Exploring the range [0-0.07] could be theoretically motivated by assuming some entropy 
production mechanism, like the decay of some other particles into the SM thermal bath. In 
that case, more ingredients are needed, and physical effects on the CMB and LSS observables 
are a bit different. Indeed, in the limit of very small DR density, one can reach a new regime 
in which the effect of the drag of DM on DR can be as relevant as the drag of DR on DM. 
We defer the study of this other class of models to a future publication. 


Parameter 

CMB+BAO 

CMB+LSS 

CMB+BAO 

CMB+BAO 




+LSS 

+LSS+77 0 

lOOwf, 

9 oq^+0.024 

z.zou_ 0 0 26 

9 91 n+0.029 

z.zi j_o.o41 

9 oon+0-021 
z.zzu-o 0 25 

9 9Q/I+0.025 

z.zo^—o 026 

Wdm 

0 1244 +0 ' 0021 
u.izuu_ 0 0040 

0 1256 +0,0034 

U.izoo_ 0 . 0047 

0 1940 + 0.0023 
u.izuy _ 0 0049 

0 1 974+0.00 40 

u.izm_ooo60 

AAfiuid 

< 0.58 

< 0.71 

< 0.67 

< 0.59 

10 7 T 0 [Mpc -1 ] 

< 1.54 

1 74+0.57 

1 (jr:+0.42 
1 .°°_0.44 

1 fjq+O-43 
i ’ uy -0.48 

Hq [km/(s Mpc)] 

69.it?| 

69-0+4 

69.ll?! 

70.2+ 3 

10 9 A s 

9 990 + 0.079 

Z.ZZU_Q 081 

9 90^+0.063 

Z.ZUO—Q Q76 

9 9nc+0.063 

Z.ZU<J_Q Qgg 

9 91 't+0.062 

Z ‘ Z1/ -0.070 

n s 

0-9709l°;°°| 

0-9762+™ 7 ? 

0 q74«+0.0051 
u.y (oo_Q 0055 

0.9796l 9 ;“l 

Treio 

0-084+°} 8 

0.078+;°i9 

0 o79 +0 ' 015 

U.U( J_ 0 0 15 

0.08278:811 


O SOSS+0-0082 

U.oUoo_o 0083 

0 qn+0.019 

u.oioU-Q.mg 

0 3097 + °- 0085 
u.ouu 1 _ 0 . 0083 

n qo'vi+o.oo^o 

U.OU9Z_ 0 0083 

O' 8 

O SI 1 +0-026 

U.511—0 019 

0 760 + °'° 17 
u. mu_ 0 .oi9 

0 7(59 +0 - 011 

u. mz^o.on 

0 766 +0011 

U. (OO^o.oil 

Ay 2 / ACDM 

0 

-9.6 

-11.4 

-12.7 


TABLE I: Mean value and 68 %CL confidence interval (or, in a few cases, 95%CL upper limit) for 
the eight parameters of our model (assuming flat priors) and two derived parameters. The last line 
shows the minimum y 2 value compared to that of ACDM with the same data. Note that 10 7 ro in 
Mpc -1 is equal to 10 21 To in s _1 to 3% accuracy. 


3 In principle, if we were using the Multinest algorithm instead of the Metropolis-Hastings algorithm, we 
could also report the Bayesian evidence ratio between the two models. This is beyond the scope of this 
paper. 
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B. Reconciling cosmological data sets 


Our results are summarized by Table [TJ The most striking facts are, first, that our model 
can reconcile CMB, BAO and LSS data, and even the Hq measurement of [16]; and second, 
that when at least CMB and LSS data are included in the fit, the minimum effective y 2 
decreases by a substantial amount when going from the ACDM model to our model: Ay 2 = 
-9.6 for CMB+LSS and Ay 2 = -11.4 for CMB+BAO+LSS. 



FIG. 4: 68% and 95% CL contours for (as, Ho) and (c78,O m ): first, for the ACDM model and 

CMB+BAO data (green); next, for our model and CMB+BAO data (black), CMB+LSS data 
(blue), CMB+BAO+LSS data (red). This figure can be compared with Fig. 33 of Planck 2015 [TJ, 
to show a clear difference between our model and all the massive active/sterile neutrino models 
used in that figure: our model can explain a lower as without requiring at the same time a lower 
Ho or a higher (on the contrary, it is compatible with higher Hq values). 

A good way to appreciate these results is to look at the (as, H 0 ) and (crg,O m ) contours 
shown in Figure|4} The CMB+BAO results for ACDM are shown in green. These results are 
notoriously in 3-4cr tension with LSS data, which require at the same time a lower as and a 
similar O m , and in 2-3er tension with the high value of Hq from [16] . The CMB+BAO results 
for our model are shown in black/grey. The comparison of the green and black contours 
makes the point. Our model is compatible with much lower values of a 8 for the same range 
of Q m values. It is also compatible with much larger H 0 values. It is worth stressing that 
this represents a crucial difference between our model and more traditional models featuring 
extra relativistic or massive relics (like sterile neutrinos) in combination with massive active 
neutrinos. These models have been invoked by Ref. [33] to reconcile tensions between CMB, 
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LSS and Ho data. The Planck 2015 paper has shown that this does not work well anymore 
with recent CMB and BAO data. Ref. [I] shows that the improvement is only of the order 
of Ay 2 ~ 1 between the ACDM model and these models. Fig. 33 of [T] provides a clear 
interpretation of this result. In models with extra massless/massive neutrinos, parameter 
correlations are such that a reduction of erg requires higher values of fl m and smaller values 
H 0 , which exacerbates tensions. We emphasise that this is not the case in our model, as 
can be seen from the black contours in Figure |4} the effect of the interaction rate allows 
for a lower erg without correlated effects on Q m , and H 0 . Hence we do obtain a substantial 
improvement in y 2 . 

The blue contours in Figure [4] correspond to CMB+LSS results for our model. A discrep¬ 
ancy between black and blue contours would reveal a tension between BAO and LSS data 
(given the CMB data). Most extensions of the minimal ACDM model studied so far lead 
to such a tension. Instead, for our model, the joint two-dimensional contours overlap at the 
l.lcr level, and the 1-dimensional posteriors for erg are compatible at the 1.4er level. This 
shows that our model actually reconciles CMB, BAO and LSS data. 

The joint contours for CMB+BAO+LSS are shown in red. Note that they are ~3er away 
from the green contours of the ACDM model, which is consistent with the typical level 
of tension between CMB+BAO and LSS data in the ACDM case, and with the fact that 
Ay 2 = —11.4 for CMB+LSS+BAO. The preferred parameter values for this data set are 
(7 g = 0.762 ^q;q}(, Q m = 0.3097^o;gog3, H 0 = 69.ll?;® km/s/Mpc, all at 68%CL. The latter 
result on Ho is compatible with the measurement of [16] at the ler level. Hence it is legitimate 
to combine CMB+BAO+LSS data with this measurement. In that case, the minimum y 2 
decreases by 12.7 with respect to ACDM. 


C. Physical interpretation of the results 

Figure [5] shows the posterior probability of the parameters of our model for the four 
combinations of data that we investigated. Since the dark matter-dark radiation interaction 
rate suppresses the matter power spectrum on small scales, this rate is compatible with zero 
for CMB+BAO data, and 3-4cr away from zero as soon as LSS data is introduced. The 
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FIG. 5: Posterior probabilities for the eight parameters forming the basis of our model and for two 
derived parameters (J) m , as), for CMB data combined with BAOs (black), LSS (blue), BAO+LSS 
(red), BAO+LSS+i/o (yellow). See the text for details on parameter definitions and priors, and 
for the precise content of each dataset. 


figure also shows that adding an H 0 prior has very little effect, excepted on H 0 itself, and 
also on AiVfluid, due to the correlation between these two parameters discussed in the next 
paragraph. 



A A fluid, 


ANfi U id 


io + 8 r 0 [ m P c - 1 ] 


FIG. 6: 68% and 95% CL contours for (AAfl u ;d, To), (AiVfl u jd, Hq), (Fo, erg), with CMB+BAO 

data (black), CMB+LSS data (blue), CMB+BAO+LSS data (red). 

The relation between the observable parameters (erg, H 0 ) and the fundamental parameters 
(r 0 , AiVfl U id) is better illustrated by figure [6] The fact that the dark matter-dark radiation 
interaction has the effect of reducing the small-scale matter power spectrum is directly 
responsible for the strong correlation between erg and Tq. Concerning the correlation between 
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H 0 and AAfl U i d , a few comments are in order. At the level of background cosmology, AAfl uid 
has a model independent effect (similar for our interacting dark radiation or for ordinary 
decoupled relativistic relics, like very light sterile neutrinos). It is well-known that an increase 
in the radiation density is compatible with CMB+BAO data provided that the characteristic 
redshifts z eq , z\ of radiation/matter and matter/A equality remain constant, and in our basis 
of parameters, this means that Hq must increase. This leads to a parameter degeneracy 
between AA fiuid and H 0 , that we clearly see in figure [6j This degeneracy is not perfect, 
due to an enhancement of the Silk damping effect when H 0 increases, and to more subtle 
perturbation effects. The latter effect is more model dependent: perturbation effects are less 
pronounced for interacting dark radiation, since it does not tree-strea m like extra massless 
relics. Nevertheless, we find bounds on AAfl u ; d for our dark radiation model which are 
comparable to those in the case of standard extra relativistic relics. For instance, with 
CMB+BAO data, we get AA fluid < 0.58 (95%CL); while with the same CMB and BAO 
data, Planck 2015 P found AA eff < 0.61 (95%CL) (with a different prior, equivalent to 
—3 < AA eff < 3, instead of our prior AA fluid > 0.07). 

Note that Planck data is more compatible with the presence of three standard free- 
streaming neutrinos than with an equivalent amount of radiation with effective parameters 
(Cg, c7 g ) 7 ^ (1/3,1/3) |H SUES]. This just means that we do observe the free-streaming effect 
of the SM active neutrinos in the CMB spectrum. But when we allow for extra species, we 
can accommodate some amount of dark radiation even if it does not free stream. 

In summary, in our fits of the data, we observe two simple and independent effects: the 
rate r 0 helps to reduce cr 8 , while the relic density of dark radiation, parametrized by AAfl uid , 
helps to increase H 0 . In the class of models studied here, which have AAfluid > 0.07, there 
is no obvious correlation between these two effects, as shown by the left plot in Figure [ 6 } 
The interaction rate can efficiently suppress a 8 even when AAfl u ; d = 0.07, the minimum 
attainable value in our particle physics models. The best fit without H 0 prior corresponds 
to the Abelian model in which the light fermions do not thermalize before the dark sector 
decouples. The best fit with the H 0 prior is close to the SU( 3) non-Abelian model, and in 
any case SU(N) with A = 4 or higher is clearly disfavored. Since there is no significant 
correlation between the effects of r 0 and AAfl uid , we get evidence for a non-zero interaction 
rate solely from LSS data, independently of the value of H 0 , which is still subject to strong 
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observational uncertainties. 


V. CONCLUSIONS 


In this paper we showed how models of dark matter interacting with dark radiation can 
reconcile CMB and LSS data. This requires a smooth reduction in the growth of modes that 
entered the horizon before matter domination, this includes the modes associated with the 
characteristic scales relevant to as- We showed that such a smooth reduction of power arises 
if the momentum transfer between DR and DM is very small and scales as T 2 , such that it 
can be equally important throughout the radiation dominated epoch but small enough to 
avoid dark acoustic oscillations. To the best of our knowledge the class of models studied in 
this paper are the first to reconcile LSS data with CMB without degrading the goodness of 
the fit to the latter. 

We have modified the Boltzmann code CLASS to include interacting DM and DR, and 
parametrized the energy density in DR by the equivalent effective number of neutrinos, 
AiVfluid, and the interaction strength by the momentum transfer today To- Motivated by 
particle physics models in which dark matter and dark radiation were in thermal equilibrium 
with the Standard Model in the early universe we have focused on a minimum value of 
AiV fluid > 0.07. 

We use the code MontePython (2BJ to fit the model to CMB, BAO and LSS data 
and found that our model improves y 2 by 11.4 compared to ACDM. The best fit has To = 
( 1 . 65 + 0 : 44 ) 10" 7 Mpc -1 (~ 1.6 x lO-^s" 1 ), Q m = 0.3097±g;gg||, H 0 = 69.1±$;| km/s/Mpc, all 
at 68 %CL, and AlVfiuid < 0.67 at 95%CL. We find 3 — 4 a evidence for non-zero momentum 
transfer rate between DM and DR. We also find that our model can accommodate larger 
values of H 0 , which have been reported by direct measurements |16j . This is a significant 
improvement over ACDM, for which an increase in H 0 exacerbates the tension with LSS 
data because it leads to an even larger prediction for erg. 

Dark matter drag also predicts a smooth suppression of the matter power spectrum at 
scales smaller than the ~ 8 Mpc relevant for <jg. Thus there may be observable effects from 
the Lyman-ct forest and perhaps also at even smaller scales where there are difficulties with 
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reconciling observations and simulations of DM halos (see [2] 3. 36] and also [37J138]). 

In the near future there will be a wealth of new data, especially for LSS (DES, LSST, 
Euclid, etc.). If the new data sets continue to favor an interacting DM and DR scenario over 
ACDM, this could be the first clear experimental evidence for non-gravitational interactions 
of DM. 
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